
********************************************************************************
$ontext
This is the numerical bottom-up optimization model used in
St�ckl, F. and A. Zerrahn (2020): Substituting Clean for Dirty Energy: A Bottom-Up Analysis

This model is derived from
The Dispatch and Investment Evaluation Tool with Endogenous Renewables (DIETER).
Version 1.3.0, November 2018.
Written by Alexander Zerrahn and Wolf-Peter Schill.

This work is licensed under the MIT License (MIT).
For more information on this license, visit http://opensource.org/licenses/mit-license.php.
Whenever you use this code, please refer to http://www.diw.de/dieter.
We are happy to receive feedback under azerrahn@diw.de and wschill@diw.de.
$offtext
********************************************************************************




*****************************************
**** Scenario file                   ****
****                                 ****
*****************************************

Parameter
m_exog_p(n,tech)
m_exog_sto_e(n,sto)
m_exog_sto_p_in(n,sto)
m_exog_sto_p_out(n,sto)
;

m_exog_p(n,tech) = technology_data(n,tech,'fixed_capacities') ;
m_exog_sto_e(n,sto) = storage_data(n,sto,'fixed_capacities_energy');
m_exog_sto_p_in(n,sto) = storage_data(n,sto,'fixed_capacities_power_in');
m_exog_sto_p_out(n,sto) = storage_data(n,sto,'fixed_capacities_power_out');


N_TECH.lo(n,tech) = 0 ;
N_TECH.up(n,tech) = 0 ;

N_TECH.up(n,'wind_on') = m_exog_p(n,'wind_on') ;
N_TECH.up(n,'wind_off') = m_exog_p(n,'wind_off') ;
N_TECH.up(n,'PV') = m_exog_p(n,'pv') ;
N_TECH.up(n,'CCGT') = m_exog_p(n,'CCGT') ;

N_STO_P_IN.up(n,sto) = 0 ;
N_STO_P_OUT.up(n,sto) = 0 ;
N_STO_E.up(n,sto) = 0 ;

N_STO_P_IN.up(n,'sto1') = m_exog_sto_p_in(n,'sto1') ;
N_STO_P_OUT.up(n,'sto1') = m_exog_sto_p_out(n,'sto1') ;
N_STO_E.up(n,'sto1') = m_exog_sto_e(n,'sto1') ;

N_STO_P_IN.up(n,'sto5') = m_exog_sto_p_in(n,'sto5') ;
N_STO_P_OUT.up(n,'sto5') = m_exog_sto_p_out(n,'sto5') ;
N_STO_E.up(n,'sto5') = m_exog_sto_e(n,'sto5') ;

N_STO_P_IN.up(n,'sto7') = m_exog_sto_p_in(n,'sto7') ;
N_STO_P_OUT.up(n,'sto7') = m_exog_sto_p_out(n,'sto7') ;
N_STO_E.up(n,'sto7') = m_exog_sto_e(n,'sto7') ;

N_TECH.up(n,'bio') = 0 ;

nu_x = 50000 ;
delta_x(n) = 0.25 * sum( h , delta(n,h) ) ;